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Abstract. We extend the formalism of dark matter directional detection to arbitrary one- 
body dark matter-nucleon interactions. The new theoretical framework generalizes the one 
currently used, which is based on 2 types of dark matter-nucleon interaction only. It includes 
14 dark matter-nucleon interaction operators, 8 isotope-dependent nuclear response functions, 
and the Radon transform of the first 2 moments of the dark matter velocity distribution. We 
calculate the recoil energy spectra at dark matter directional detectors made of CF4, GS2 and 
^He for the 14 dark matter-nucleon interactions, using nuclear response functions recently 
obtained through numerical nuclear structure calculations. We highlight the new features 
of the proposed theoretical framework, and present our results for a spherical dark matter 
halo and for a stream of dark matter particles. This study lays the foundations for model 
independent analyses of dark matter directional detection experiments. 
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1 Introduction 

The Earth’s motion with respect to the galactic rest frame is expected to produce a flux 
of dark matter particles across the surface of the planet [1]. Fixed target experiments can 
in principle detect these particles, if they scatter and deposit energy in low-background 
detectors [2-4]. 

The angular distribution of nuclear recoil events from dark matter-nucleus scattering at 
fixed target experiments is not isotropic, as the Earth’s motion selects a preferred direction in 
the sphere of recoil directions [5] . Recoil events are mainly expected in the direction opposite 
to the observer’s motion, or in a ring around it for large dark matter particle mass to recoil 
energy ratios [6]. 

Dark matter directional detectors are designed to measure anisotropies in the distribu¬ 
tion of nuclear recoil events at fixed target experiments. Specifically, their task is to measure 
the recoil momentum vector of nuclei scattered off by dark matter particles from the local 
galactic population [7] . The sense of the recoil momentum vector is also called the “head-tail” 
of the recoil track. 

Dark matter directional detection experiments currently in a research and development 
stage are DRIFT [8, 9], MIMAC [10, 11], DMTPC [12, 13], NEWAGE [14, 15] and D3 [16]. 
They adopt diffuse gas detectors and time projection chambers to reconstruct the nuclear 
recoil tracks. Alternative approaches include the use of nuclear emulsions [17], dark matter- 
electron scattering in crystals [18] and DNA detectors [19]. In the near future, directional 
detection experiments should reach the sensitivity to detect standard weakly interacting 
massive particles (WIMPs), as those predicted by SUSY [20-23] or extra-dimensional theo¬ 
ries [24-27]. 
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The recoil energy spectrum expected at dark matter directional detectors has so far 
been computed for two types of dark matter-nucleon interaction only, e.g. [28-37]. One is 
the familiar spin-independent interaction. The other one is the well-known spin-dependent 
dark matter-nucleon interaction. At the quantum mechanical level, both interactions are 
independent of the momentum transfer operator, and of the dark matter-nucleon relative 
velocity operator. 

Though assuming momentum and velocity independent dark matter-nucleon interac¬ 
tions is a reasonable first approximation, it only provides a limited, and to some extent biased 
description of the actual complexity of the dark matter-nucleus scattering [38]. The most 
general non-relativistic effective theory for one-body dark matter-nucleon interactions pre¬ 
dicts 14 Galilean invariant dark matter-nucleon interaction operators, each with an isoscalar 
and an isovector component [39, 40]. On purely observational grounds, there is no reason to 
neglect any of the 14 interaction operators, which motivates the exploration of more general 
approaches. 

In this work we extend the formalism of dark matter directional detection to arbitrary 
one-body dark matter-nucleon interactions. The ultimate goal of this study is to set the bases 
for fully model-independent analyses of future dark matter directional detection experiments. 
Non-relativistic dark matter-nucleon interaction operators were previously explored in the 
context of dark matter direct detection in [38-59] , and studying the dark matter capture by 
the Sun in [60-67]. 

The paper is organized as follows. In Sec. 2 we introduce the non-relativistic effective 
theory of dark matter directional detection. The theory is carefully explored in Sec. 3 for 
three benchmark detectors and two astrophysical configurations. We summarize our results 
in Sec. 4, and list useful equations in the Appendixes. 

2 Effective theory of dark matter directional detection 

In this section we extend the formalism of dark matter directional detection to arbitrary 
one-body dark matter-nucleon interactions. 

2.1 Dark matter-nucleus scattering in effective theories 

We start with a brief review of the effective theory of one-body dark matter-nucleon inter¬ 
actions [39]. It provides a general framework for the study of dark matter scattering from 
target nuclei at directional detectors. 

Under the assumption of one-body dark matter-nucleon interactions, the Hamiltonian 
density 

= E E E *f.) (2-1) 

i=l r=0,l k 

describes the most general short-range dark matter-nucleus interaction [39]. In Eq. (2.1), the 
sum over i = I,..., A, where A is the target nucleus mass number, reflects the assumption 
of one-body interactions. The 14 dark matter-nucleon interaction operators depend on 
the momentum transfer operator q, on the relative transverse velocity operator v-*-, and on 
the dark matter particle and nucleon spin operators, and Sw, respectively [39]. We list 
them in Tab. 1. Introducing the 2x2 matrices in isospin space = 1 and = T 3 , where 
r 3 is the third Pauli matrix, we relate the isoscalar and isovector coupling constants, 
and respectively, to the coupling constants for protons (c^) and nucleons (c^) as follows: 
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^3 = iSN ■ X V^) 

O 4 = • Sjv 

^5 = iS^ ■ X v^) 

= (Sx • ^) (Sjv • ^) 
dr = Sn-v^ 
ds = S^- 


a 9 = iS^- S^x-^ 


mjv 


Oio = ^Sat • ^ 
On = ^ 

di2 = s 

Oi3 = ii Sx 



Ol4 = ^(S,•^ 

Oi 5 = -fs^-;;^ 


q 

mjv 


Table 1. Non-relativistic quantum mechanical operators defining the general effective theory of one- 
body dark matter-nucleon interactions. All operators have the same mass dimension, and m^v is the 
nucleon mass. Here and in the next sections we omit the nucleon index (i). 


= (c° -|- c\)l2, and = (c° — c],)/2. Defined in this way, and c\ have dimension mass 
to the power of —2. The Hamiltonian density (2.1) admits the following coordinate space 
representation 


{ A A 

^ S(r - Fj) -b ^ • di 6{r - Vi 

i=l 


2 = 1 


A 


+ 


2mjv 

+ y^E- — 

^ A 9 


i=l 

A 


i=l 


2m jq 


' i^rS{r — Ti) — i5{r — rj)^r 

X ai6{r - Ti) + (5(r - rj) Ui x 




( 2 . 2 ) 


It depends on the Pauli matrices that represent the zth-nucleon spin operator a,, on the dark 
matter-nucleus relative distance r, and on the fth-nucleon position in the nucleus center of 
mass frame r* [39]. 

Eq. (2.2) describes the dark matter coupling to the nuclear vector charge (first term 
in the first line), to the nuclear spin-current (second term in the first line), to the nuclear 
convection current (second line), and hnally, to the nuclear spin-velocity current (last line). 
The exact nature of these couplings is encoded in the four operators Iq, I 5 , and They 
are defined as follows 


Iq = c[ + i 


’ X y) • Sx C 5 + • Sx 4 + • Sx Cn 

V rriM J mN 


1 


X v;^ C3 -b C4 -b —-^ -8x06+ Vt C7 + x Cg -b c[o 

mjv ^ miq mN ^ ruN ^ tun 


+ X S, • s, C [3 + . S, Cl 4 + ^ X V# ^ . S, C 35 

\lr = t^_xS^cl-S^cl 


mN 


Hricr .qq^r 

— C3 + C12 - — X C33 - i—— ■ Sx Ci5 
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(2.3) 


The four operators Iq, I 5 , V^, and depend on the dark matter particle spin operator 
S^, on the momentum transfer operator q, and on the operator = v"*- — v]^, where 
= i^rS{r — rj) — i6{r — rj)^r- Analogous coordinate space representations for q 
and v;); can be found in [ 68 ] . 

The Hamiltonian density in Eq. (2.2) leads to the following transition probability for 
dark matter-nucleus scattering 


(I Ad t \ ) spins — 


dvr 


2J+1 
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+ 


m 
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E 


Rl 




y¥\^jwr(y) 


(2.4) 


where Mt is the transition amplitude, J is the nuclear spin, and the angle brackets denote 
a sum (average) over the final (initial) spin-configurations. The index k in Eq. (2.4) extends 
over the 8 nuclear response functions [y) defined below in Sec. 2.2 together with the 
new kinematical variable y. 

We refer to the 8 functions in Eq. (2.4) as dark matter response functions. They 
are quadratic in matrix elements of the operators in Eq. (2.3), and depend on and 

on where v is the dark matter-nucleus relative velocity, is the 

nucleon mass, and yx the reduced dark matter-nucleus mass. We list the functions in 
Appendix A. 

The differential cross-section for dark matter scattering from nuclei of type T and mass 
niT can finally be written as follows 


do-r(g^,'c^) 

d (/2 


1 


(jAd'T’j )spins ) 


(2.5) 


which for arbitrary interactions is a function of the momentum transfer, and of the dark 
matter-nucleus relative velocity. 


2.2 Nuclear response functions and target materials 

The 8 nuclear response functions in Eq. (2.4) arise from a multipole expansion of the nuclear 
charge and currents in Eq. (2.2), and are defined as follows 

Wl];{y) = Y,{J^T,MT\\ aL-ri<l) \\J,T, Mt){J,T, Mt\\ hy{q) UT,Mt) . (2.6) 

L 

The multiple expansion index L must be less then 2J. In Eq. (2.6), we label the state 
I J, T, Mt) using the nuclear spin J, the nuclear isospin T, and the isospin magnetic quantum 
number Mt- The functions are therefore quadratic in nuclear matrix elements reduced 
in the nuclear spin magnetic quantum number. Assuming the harmonic oscillator basis for 
single-nucleon states, the nuclear response functions depend on y = {bq/2)‘^ only, where 

b = y^41.467/(45A“^/^ — 2f)A~‘^l^) fm . 
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(2.7) 








The nuclear response operators aLM-,T and hr Eq. (2.6), can each be one of the following 

operators 
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( 2 . 8 ) 


where Miuiqri) = jLiq^YLuAn), and M^iq^i) = jAq^YAiArt)■ The vector spheri¬ 
cal harmonics, Y^^( 0 rJ, are defined as follows 


yA'iAA = Y,{L'mlX\L'lLM)YL>^{nr^)ex , (2.9) 

mX 

where (L'mlA|L'lLM) are Clebsch-Gordan coefficients, and is a spherical unit vector 
basis. In Eq. (2.4), the index k extends over pairs a/3 of nuclear response operators, e.g. 
k = AS'. We however use the notation {y) = (y) for a = ^. 

Dark matter directional detection experiments currently in a research and development 
stage mainly exploit target materials made of CE 4 , CS 2 and ^He, or of mixtures of them. 
Eor instance, experiments like DMTPC and NEWAGE adopt detectors composed of GE 4 , 
whereas DRIET and MIMAC use, respectively, GS 2 (or a GS 2 -CE 4 mixture) and ^He (or 
CE 4 ) as target materials. 

Eor the elements ^He, and we use the nuclear response functions derived in [ 68 ] 
through numerical nuclear structure calculations. Eor ^®E we adopt the nuclear response 
functions found in [39]. Eor reference, the nuclear response functions relevant for dark matter 
directional detection are listed in Appendix B. 


2.3 Recoil energy spectra 

Eor a given dark matter-nucleon interaction operator in Tab. 1, and a given target material, 
the double differential nuclear recoil energy spectrum per unit time and per unit detector 
mass expected at a directional detection experiment is given by 


d‘^n 
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E 
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( 2 . 10 ) 
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where Er = j{2mT) is the nuclear recoil energy, w the nuclear recoil direction (|wp = 1), 
the dark matter particle mass, the mass fraction of the nucleus T in the target 
material, and py the local dark matter density. In Eq. (2.10), wt = <?/2/rT is the isotope- 
dependent minimum velocity required to transfer a momentum q from the target nucleus to 
the dark matter particle, and Ve(t) is the time-dependent Earth’s velocity in the galactic 
rest frame. The angle 0 is measured with respect to the reference direction Ve(t). Here we 
assume azimuthal symmetry of (^TZ/dE^dQ. around the direction Ve(t), i.e. dfl = 27rdcos0. 
The integral in Eq. (2.10) corresponds to the Radon transform /m of the local dark matter 
velocity distribution / in the galactic rest frame boosted to the detector frame times the 
transition probability (2.4): 

fMiwT,w) = j d{v-W-WT)f{v + Ve{t)){\MT\‘^)spmsd^V. (2.11) 


Our assumptions regarding /(v -|- Ve(t)) and Ve(t) will be discussed in the next subsection, 
where we also provide analytic expressions for the Radon transform /m, given the transition 
probability (2.4). 

In the figures of Sec. 3 we plot the differential number of recoil events around cos 9 per 
unit time and per unit detector mass: 


d7^ ^ f d^TZ 

-— 2i7r I - 

dcos6l JER>E,i, dERdn 


( 2 . 12 ) 


where Eth is the detector energy threshold. 


2.4 Radon transforms 

For /(v -|- Ve(t)) we assume a normalized Maxwell-Boltzmann distribution truncated at the 
escape velocity Uesc^ 


/(v + Ve(t)) = 


1 


iVesc(27ru2)3/2 V 2^2 

where the normalization constant A^esc is given by 

2 Uesc 


V Ve 


0 (Vesc - |v Ve 


(2.13) 


Nesc = erf 


V2, 


(7 7 ; 


TT (T')) 


exp 


2 ct 2 


(2.14) 


and is the velocity dispersion. In Eqs. (2.13) and (2.14), Uesc and Ve(t) are measured in 
the galactic rest frame, whereas v is the dark matter-nucleus relative velocity in the detector 
frame. 

Eq. (2.13) is a standard assumption for the local dark matter velocity distribution. 
Though a first approximation only, it allows to analytically compute /M(w'r,w) for a given 
transition probability (|AlTp)spins- We leave the exploration of velocity distributions self- 
consistently generated by Eddington’s inversion method [69], or by its anisotropic exten¬ 
sions [70, 71], for future work. 

Inspection of Eq. (A.l) shows that for arbitrary one-body dark matter-nucleon inter¬ 
actions, two independent Radon transforms appear in Eqs. (2.10) and (2.12). One, for 
(|A4rP)spins oc the other one , for (|A4rP)spms oc v'^. 
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Within the approximation (2.13), we can analytically calculate both and We 

find 

j (5(v-w-'u;r)/(v + ve(t))d^v 

1 J (tCT + |Ve(f)| cos 0)^ 

" Aresc(27rcT2)l/2 [ 2^2 



-exp , (2.15) 


/m ('»^T,w) = J (5(v. w-'u;r)/(v + Ve(t))u^d^v 

(ref. + |ve(t)|^ sin^ 0 + 2cr 
cos 6 {wt + I Ve (t) I cos 9) 

(2.16) 

for {wt + |ve(t)| cos 9) < Uesc, and zero otherwise. We have verihed the validity of Eqs. (2.15) 
and (2.16) using the Radon transform inversion theorem, which under minimal regularity 
conditions relates a function g to its Radon transform g as follows 

5(v) = -^ /5"(v. w,w)dnu;. (2.17) 

In Eq. (2.17), g''{x, w) = d‘^g{x, w)/9x^, x = v-w, and the integral is over all recoil directions 
w. Using Eq. (2.17), we were able to re-derive /(v+Ve(t)) and u^/(v+Ve(t)) from Eq. (2.15) 
for f^\wT,v/) and Eq. (2.16) for f^\wT,y^), respectively. Our Eqs. (2.15) and (2.16) allow 
to evaluate Eqs. (2.10) and (2.12) for all one-body dark matter-nucleon interaction operators 
in Tab. 1. 

The Earth’s velocity Ve(f) in the equations of Secs. 2.3 and 2.4 is the vectorial sum of the 
local standard of rest velocity, of the solar motion, and of the Earth’s velocity with respect 
to the Sun. Here we consider a local standard of rest velocity of 220 km s“^, a solar motion 
of (11.1,12.2, 7.2) km s“^, and neglect the remaining (time-dependent) contribution to Ve(t). 
As for the Maxwell-Boltzmann approximation (2.13), this is a standard assumption regarding 
Ve(t), but it is not free from uncertainties. We refer to [72-75] for an exhaustive discussion 
on this subject. If not otherwise specified, we assume = 156 km s“^, p-^ = 0.4 GeV cm“^, 
and Uesc = 533 km s“^ [73, 76, 77]. 

3 Phenomenology 

In this section we study the phenomenology of the effective theory of dark matter directional 
detection defined in Sec. 2. The main observable we are interested in is the differential number 
of recoil events around cos0 per unit time and per unit detector mass, namely dR/dcos0 
in Eq. (2.12). Using the Radon transforms of Sec. 2.4, and the nuclear response functions 
dehned in Sec. 2.2, and given in Appendix B, we now compute dR/dcos0 for all dark matter- 
nucleon interaction operators in Tab. 1. In Secs. 3.1 and 3.2 we focus on scattering events 
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Figure 1. Differential number of recoil events around cos0 per unit time and per unit detector mass 
as a function of cosd, i.e. d_R/dcos0 in Eq. (2.12). We report di?/dcos0 for directional detectors with 
head-tail discrimination made of CS 2 , CF 4 and ^He. Solid, dashed and dot-dashed lines correspond 
to equal to 10 GeV, 100 GeV and 1000 GeV, respectively, as shown in the legends. The top 
panels refer to the coupling constants and c}, whereas the bottom panels correspond to C 4 and C 4 , 
respectively. Here we assume a truncated Maxwell-Boltzmann distribution for the Milky Way dark 
matter particles. 


induced by dark matter particles from the Milky Way dark matter halo. As an example of 
an alternative astrophysical configuration, in Sec. 3.3 we consider a hypothetical stream of 
dark matter particles moving in the direction of — Ve(t). If not otherwise specified, we assume 
Eth = 0. Only in Sec. 3.4 we relax this assumption, and quantify the impact of Et^ / 0 on 
di?/dcos 9. 

In the following, the angle 9 is measured with respect to the reference direction Ve(t). 
To simplify the comparison of different figures, in all plots we assume a single coupling 
constant different from zero and equal to 10 “^/m^ at the time, where = 246.2 GeV is 
the electroweak scale. Results are presented for three benchmark dark matter directional 
experiments with head-tail discrimination made of CS 2 , CF 4 , and ^He, respectively. For 
definitiveness we consider spin- 1/2 dark matter. 
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Figure 2. Same as Fig. 1, but now for the operators and O^. 


3.1 Momentum/velocity independent operators 

We start with an analysis of the only momentum transfer and relative velocity independent 
operators, namely 0\ and O 4 . In Eq. (2.2), the operator 0\ contributes to the nuclear vector 
charge through the nuclear response operator which in the low momentum transfer 

limit measures the nuclear mass number. The operator O 4 contributes to the nuclear spin 
current through and which in the low momentum transfer limit measure the 

nucleon spin content of the nucleus. 

The cross-sections for dark matter-proton and dark matter-neutron scattering derived 
from the operator Oi are 

,,2 UO I U|2 ,,2 |„0_U|2 

.^SI Tl+Cil SI fijY |Ci Cil 

= EW-z-’ = EW-z-’ 


respectively. Similarly, from the operator O 4 one obtains the cross-sections for dark matter- 
proton and dark matter-neutron scattering 



^l/JxOx + 1) |c4 + ciP 
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SD 


+ 1) |C4 - cjp 
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Figure 3. Same as Fig. 1, but now for the operators Oq and Oj. 


respectively. Here hn = /{^x + '^n) is the reduced dark matter-nucleon mass, and 

= 1/2 is the dark matter particle spin. 

Eqs. (3.1) and (3.2) connect the formalism reviewed and developed here, to standard 
analyses of dark matter directional detection, where = c\ and C 4 = c\ is commonly 
assumed. In contrast, here we consider c)*, c}, C 4 , and c\ as independent parameters, and 
compute di?/dcos 0 accordingly. 

Fig. 1 shows the differential rate di?/dcos0 as a function of cos0 for the operators Oi 
and O 4 . In the left panels we focus on the isoscalar coupling constants c? and C 4 , whereas 
the right panels are devoted to the isovector coupling constants c\ and c\, respectively. In 
the plots different colors and lines refer to the target materials, and to the values of the dark 
matter particle mass reported in the legends. 

The isoscalar component of the operator Oi can be probed by detectors made of CS 2 , 
CF 4 , and ^He. In contrast, a detector composed of CS 2 would be insensitive to the isovector 
component of the operator Oi. 

Regarding the operator O 4 , we find that this interaction cannot be detected using CS 2 
as a target material, as both and have spin zero. We also observe that the isoscalar 
and isovector components of the operator O 4 generate a similar differential rate dii/dcos 0 . 
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Figure 4. Same as Fig. 1, but now for the operators Os and Og. 


as the nuclear response functions differ by about a factor of 2 for 

and ^He (see appendix B). 

3.2 Momentum/velocity dependent operators 

Next, we focus on the dark matter-nucleon interactions in Tab. 1 depending on the momentum 
transfer operator q, on the dark matter-nucleus transverse relative velocity operator v-*-, or 
on both. 

Contrary to Oi and O 4 , for these operators the differential rate di?/dcos0 can depend 
on nuclear response functions that do not appear in the theory of electroweak scattering 
from nuclei, and are specific to one-body dark matter-nucleon interactions. An example of 
nuclear response functions that are not generated by the operators Oi and O 4 is W^T/ , which 
is different from zero for and The corresponding nuclear response operator 

measures the content of nucleon spin-orbit coupling in the nucleus. Therefore, 
W^T/ can be large for isotopes with non fully occupied single-nucleon orbits of large angular 
momentum [39]. 

Figs. 2, 3, 4, 5, 6 and 7 show the differential rate di?/dcos0 as a function of cosO for the 
interaction operators in Tab. 1 that depend on q and/or v-*-. Different colors correspond to 
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Figure 5. Same as Fig. 1, but now for the operators Oio and On. 

distinct target materials. As in the previous section, we compute di?/dcos0 for directional 
detectors with head-tail discrimination made of CS 2 , CF 4 and ^He. Solid, dashed and dot- 
dashed lines denote a dark matter particle mass of 10 GeV, 100 GeV and 1000 GeV, 
respectively. 

The amplitude of the differential rate di?/dcos0 for a given interaction operator Ok 
reflects the number of q and v-*- operators that multiply in Eq. (2.3). However, it also 
crucially depends on the nuclear response functions in Appendix B integrated over energies 
larger than Eth- 

As an example of the importance that nuclear response functions can have in the calcu¬ 
lation of di?/dcos0, let us focus on the operators O 4 = ■ Sn and On = iS^ ■ q/mw- The 

former generates a nuclear spin current, the latter a nuclear vector charge. We find that for 
GF 4 detectors, the isoscalar component of the operator On contributes to di?/dcos0 with 
the same strength of the isoscalar component of the operator O 4 , though O 4 is momentum 
independent. For ^^G and F, the response function is larger than and W^,, and it 
compensates for the momentum suppression characterizing the operator On. Interestingly, 
in the literature the operator On is much less explored than the familiar spin-dependent 
operator O 4 . 
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Figure 6. Same as Fig. 1, but now for the operators O 12 and O13. 


In Figs. 2, 3, 4, 5, 6 and 7, we also observe that (besides 2 exceptions discussed below) 
the spin-dependent operators in Tab. 1, i.e. Stv dependent, do not produce any signal at 
CS 2 detectors, as both and have spin 0. The operators O 12 = • (Sat x v-*-) 

and Oi 5 = —(S^ • q/mAr)[(SAr x v-*-) • q/ruAr] constitute 2 interesting exceptions. Though 
Sat dependent, they do not only induce a dark matter coupling to the nuclear spin. They 
also generate the nuclear response function , which is different from zero for and 
Interestingly, detectors made of CS 2 cannot probe any dark matter-nucleon isovector 
coupling. 

We conclude this section with a brief discussion on the dependence of the differential 
rate (2.12). The rate di?/dcos0 depends on through the minimum velocity wt, and the 
dark matter number density {^py^jra^. The factor [pylrriy) in Eq. (2.12) is the same for all 
operators in Tab. 1. The dependence on rriy of di2/dcos 9 is therefore non trivially determined 
by the operator-dependent double integral in Eq. (2.12). Eor a given interaction operator, 
the dependence of di?/dcos0 on rriy can be qualitatively understood from a plot of the total 
scattering rate as function of rriy. Fig. 8 shows the total rate i?, defined as 


R = 


f+i 


dcos 9 


/-I 


dR 

dcos 9 ’ 


(3.3) 
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Figure 7. Same as Fig. 1, but now for the operators O 14 and O 15 . 


as a function of rriy. for all operators in Tab. 1. We find that for operators with differential 
cross-sections (2.5) independent of q, or depending on q trough only, the total rate R 
grows with up to ~ 10 GeV, and then it decreases for larger masses. This applies 
to the operators Oi, O 4 , O 7 , and Os- For operators with differential cross-sections explicitly 
depending on momentum, i.e. not through vj^ only, R grows up to of the order of 100 
GeV, and then it decreases for larger values of m^. The dependence observed in Fig. 8 
reflects in Figs. 2, 3, 4, 5, 6 and 7 accordingly. 

3.3 Changing astrophysical assumptions 

In the limit (T^/| Ve| —)• 0, the equations derived in Sec. 3.3 can also be applied to the directional 
detection of dark matter particles from galactic streams. Here we focus on a hypothetical 
stream of dark matter particles moving at a speed of 300 km s“^ in the direction of — Ve, 
with velocity dispersion = 5 km s“^ and py- = 0.3 GeV cm“^. 

Fig. 9 shows the differential rate di?/dcos0 as a function of cos 9 for the hypothetical 
dark matter stream described above, and for all dark matter-nucleon interaction operators 
in Tab. 1. Here we consider CF 4 as a target material, = 100 GeV as a benchmark dark 
matter particle mass, and Sth = 0 - 
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Figure 8 . Total rate i? as a function of for all operators in the legends. The left panel refers 
to the isoscalar couplings, whereas the right panel to the isovector couplings. For operators with 
differential cross-sections (2.5) independent of q, or depending on q trough only, the total rate 
R grows with up to ~ 10 GeV, and then it decreases for larger masses. For operators with 
differential cross-sections explicitly depending on q, R grows up to of the order of 100 GeV, and 
then it decreases for larger values of m^. 

For interaction operators independent of v*-, the differential rate di?/dcos0 has a max¬ 
imum at COS0 = —1, and it rapidly goes to zero around cos 9 = 0, as expected for the stream 
considered here. 

For interaction operators that do depend on v^, dR/dcos9 is linear in and 

it therefore exhibits a different behavior. For these operators di?/dcos0 has a maximum 

at larger values of cos 9, because of the {wt,'w) contribution to the differential rate. 

This interesting effect is not present when only the interaction operators Oi and O 4 are 

considered. The shift in the peak of di?/dcos0 is pronounced for the interaction operators 

O 5 , O 7 , Os, Ois, and £> 14 . For the v-*- dependent operators O 3 , O 12 and £>15 the effect 

2 ) 

is instead negligible, as for these operators the contribution of {wt,'n) to dR/dcos9 is 
sub-leading (see also Eq. (A.l) in Appendix A). 

3.4 Threshold effects 

We now relax the assumption Eth = 0, and study the impact of a finite experimental energy 
threshold on the differential rate di?/dcos0. Here we consider Eth = 20 keV, though lower 
energy thresholds are expected for the future [7]. We restrict the analysis to a directional 
detector made of CF 4 , and consider 10 GeV and 100 GeV as benchmark dark matter particle 
masses. 

Fig. 10 shows the ratio of the differential rate di?/dcos0 for = 20 keV to the same 
rate with Eth = 0 for all interaction operators in Tab. 1. Colors and lines refer to the 
operators reported in the legends. 

Regarding threshold effects, we find that the operators in Tab. 1 divide into 4 inde¬ 
pendent subsets. A first subset includes the operators Oi, O 4 , O 7 , and Og- These are the 
operators most sensitive to a finite energy threshold. They contribute to the differential 
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Figure 9. Differential number of recoil events around cos 9 per unit time and per unit detector 
mass as a function of cos 6 for a stream of dark matter particles moving at a speed of 300 km s“^ 
in the direction of —Vg, with velocity dispersion cr„ = 5 km s“^, and = 0.3 GeV cm“^. We 
report di?/dcos0 for a directional detector with head-tail discrimination made of CF 4 , and assume 
= 100 GeV. The left panel and the right panel refer, respectively, to the isoscalar and isovector 
components of the operators in Tab. 1, as shown in the legends. 


rate ( 2 . 12 ) through dark matter response functions either constant, or proportional to 
The remaining 3 snbsets of operators are increasingly less sensitive to threshold effects. The 
second snbset of operators involves O 5 , Og, Oiq, Oh, O 12 , O 13 , and O 14 . The dark matter 
response functions of these operators are either proportional to q^, or to (v^ contribu¬ 

tions are sub-leading). The third and fourth subsets consist of O 3 and Oq, and of O 15 alone, 
respectively. The former is characterized by dark matter response functions proportional to 
q'^ {q^ contributions are sub-leading). The latter by response functions proportional to q^ 
{q^ contributions are sub-leading). 

The dependence on the momentum transfer specific to operators in different subsets 
reflects in the 4 distinct behaviors that we observe in the top panels {rriy- = 100 GeV) and 
in the bottom panels {m-^ = 10 GeV) of Fig. 10. As expected, for low dark matter particle 
masses, and for cos 6 > — 1 , threshold effects are more pronounced. 

4 Conclusions 

We have extended the formalism of dark matter directional detection to arbitrary one-body 
dark matter-nucleon interactions. The resulting theoretical framework predicts 28 indepen¬ 
dent recoil energy spectra characterized by 8 isotope-dependent nuclear response functions, 
and by the Radon transforms of the first and second moment of the local dark matter veloc¬ 
ity distribution. In this stndy, we have analytically computed the Radon transform of the 
first 2 moments of a truncated Maxwell-Boltzmann distribution, and made use of nuclear 
response functions recently obtained in the literature through numerical nuclear structure 
calculations. 

We have computed the rate di?/dcos0 as a function of the nnclear recoil angle for the 
28 isoscalar and isovector dark matter-nucleon interactions of the theory. We have presented 
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Figure 10. Ratio of the differential rate di?/dcosd for Eth = 20 keV to the same rate with Eth = 0, 
for all interaction operators in Tab. 1. Colors and lines refer to the operators reported in the legends. 
The top panels correspond to = 100 GeV, whereas the bottom panels refer to = 10 GeV. 

our results for 3 hypothetical dark matter directional detectors with head-tail discrimination 
composed of CS2, CF4 and ^He, respectively. The relevant nuclear response functions for the 
isotopes ^He, and are listed in Appendix B. We have separately considered the 

case of dark matter particles from the Milky Way dark matter halo, and of a dark matter 
stream. 

We have found that only CF4 and ^He detectors can probe all interactions explored in 
this work (^He detectors with a typically smaller rate). CS 2 detectors are insensitive to dark 
matter-nucleon interaction operators that couple to the nuclear spin-current only, as both 
and have zero spin. Interestingly, we have found that the spin-dependent operators 
£>12 = • (Stv X V-*-) and O 15 = —(S;^. • q/mAr)[(SAr x v-*-) • q/ruAr] remain accessible to 

CS 2 detectors, in that the operators O 12 and O 15 generate a non zero nuclear spin-velocity 
current. 
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We have also found that for CF4 detectors the operator On = can contribute 

to di?/dcos 0 with the same strength of the more popular operator O 4 = S;^ • Sat, though 
O 4 is momentum independent. This intriguing result highlights the importance of nuclear 
response functions, and numerical nuclear structure calculations in dark matter directional 
detection. 

We have characterized the dependence on and on the experimental energy thresh¬ 
old ifth of the rate di?/dcos0 for all one-body dark matter-nucleon interactions in Tab. 1. 
For a given dark matter-nucleon interaction, both properties reflect the momentum transfer 
dependence of the differential scattering cross-section. We have found that the stronger the 
differential cross-section depends on the momentum transfer the less a specific interaction is 
sensitive to Flth- 

Finally, we have found that the rate dii/dcos0 for the operators O 5 , O 7 , O 13 , and 

Oi 4 has not a maximum in the direction opposite to the observer’s motion. The shift in the 
peak of dii/dcos0 is pronounced for dark matter streams. This effect is present whenever 
the second moment of the dark matter velocity distribution is quantitatively important, 
as for the operators listed above. This feature is similar to the rings found in [ 6 ] studying 
d^R/dcos 0 dEji, but it characterizes the rate di 2 /dcos 0 , which is integrated above the energy 
threshold Flth- 

Our study shows that the phenomenology of dark matter directional detection can be 
more complex and rich than in current analyses based on the familiar operators Oi and O 4 
only. The equations reported in this work can be used for model independent analyses of 
dark matter signals at directional detection experiments. 
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A Dark matter response functions 

Below, we list the dark matter response functions that appear in Eq. (2.5). The notation is 
the same used in the body of the paper. 
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B Nuclear response functions 

Below, we list the nuclear response functions relevant for dark matter directional detection 
and different from zero. For ^He, ^‘^C, and ^ 23 , we use the nuclear response functions that 
we obtained in [68]. For ^®F we adopt the nuclear response functions found in [39]. 

Helium (^He) 


W^{y) = 0.358099e"22' 
Wl}{y) = 0.0397887e"2^ 
Wj^iy) = 0.119366e-23/ 
W^{y) = 0.1193666“^*' 

Carbon (^^C) 


W^?,{y) = 0.0397887e"2^ 
Wl},{y) = 0.0397887e"2^ 
Wl?,{y) = - 0 . 03978876 - 2 ?^ 
W^},{y) = -0.03978876-2?? 


Wl9{y) = 0.07957756-2?? 
W^{y) = 0.07957756-2?? 
W^{y) = -0.07957756-2?? 
W^{y) = -0.07957756-2?? 


W^iy) = 0.5658826-2??(2.25 - yf 
Wl^,{y) = 0.04808056-2?? 

W^^„{y) = 6-2??(-0.371134 + 0.164948y) (B.2) 

Fluorine (^®F) 


W^{y) = 6 - 2 ?? (0.0662231^“^ - 1.23196y3 + 7.68018y2 _ 18.1437?/ + 14.3637) 

W}^{y) = 6 - 2 ?? (0.00683961/ - 0.043991/ + 0.103729/ _ 0.106103?/ + 0.0397887) 
W}^{y) = 6 - 2 ?? (-0.0212824/ + 0.266402/ - 1.00139/ + 1.48545?/ - 0.755986) 
W'^iy) = 6 - 2 ?? (-0.0212824/ + 0.266402/ - 1.00139/ + 1.48545?/ - 0.755986) 
W^,{y) = 6 - 2 ?? (0.00679021/ - 0.0400138/ + 0.0896114/ - 0.0903448?/ + 0.0346155) 
W^},{y) = 6 - 2 ?? (0.00716698/ - 0.0431855/ + 0.0977332y2 _ 0.0984534?/ + 0.0372496) 
VF^O(y) = e-2?? (0.00697605/ - 0.041572/ + 0.0935814?/2 - 0.0943139?/ + 0.0359084) 
W^,{y) = 6 - 2 ?? (0.00697605/ - 0.041572/ + 0.0935814?/2 _ 0.0943139y + 0.0359084) 
W^{y) = 6 - 2 ?? (0.0149629/ - 0.0867403/ + 0.19008/ - 0.186579?/ + 0.069231) 
W^{y) = 6 - 2 ?? (0.0140042/ - 0.0865149/ + 0.198218/ - 0.199543?/ + 0.0744992) 
W^{y) = 6 - 2 ?? (0.0144756/ - 0.0866714/ + 0.194127/ - 0.192953?/ + 0.0718168) 
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W^{y) = (0.0144756/ - 0.0866714/ + 0.194127/ - 0.192953?/ + 0.0718168) 

Wl^,{y) = (0.00314736/ - 0.0157368?/ + 0.019671) 

Wl},{y) = e-'^y (0.000600385/ - 0.00300193?/ + 0.00375241) 

Wl^,{y) = e-‘^y (-0.00137464y2 + 0.00687319?/ - 0.00859149) 

H/Oi (y) = e-‘^y (-0.00137464y2 + 0.00687319?/ - 0.00859149) 

W^{y) = 3.60124 x 10“®e“2?/ (2.5 - yf 
W^{y) = 0.000546075e“2^(2.5 - yf 
VF/(y) = 0.0000443458e-2y(2.5 - yf 
W^{y) = 0.0000443458e"2^(2.5 - yf 
W^^4y) = e-‘^y (0.014437/ - 0.17038/ + 0.54834y - 0.531554) 

W^^„(y) = e-'^y (0.00202642/ - 0.0115829/ + 0.0211796y - 0.012219) 

W^^„{y) = e-‘^y (-0.00463969/ + 0.02652^^ _ 0.0484926y + 0.0279765) 

W^^„{y) = e-‘^y (-0.0063055/ + 0.0744149/ - 0.239492y + 0.232161) (B.3) 

Sulfur (32S) 

W^{y) = 0.580305e"2^(5.92494 - 5.43118y + y'^f 
Wi?,{y) = 0.0765941e-2y(2.5 - yf 

W^^4y) = e-‘^y{-3.12284: + 4.11173y - 1.6721/ + 0.210827/) (B.4) 
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